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ABSTRACT 

Many studies of the Loop I magnetic superbubble place the Sun at the edges 
of the bubble. One recent study models the polarized radio continuum of Loop I 
as two magnetic shells with the Sun embedded in the rim of the 'SI' shell. If the 
Sun is in such a shell, it should be apparent in both the local interstellar magnetic 
field and the distribution of nearby interstellar material. The properties of these 
subshells are compared to the interstellar magnetic field (ISMF) and the distribu- 
tion of interstellar Fe"*" and Ca"*" within ~ 55 pc of the Sun. Although the results 
are not conclusive, the ISMF direction obtained from polarized stars within ~ 30 
pc is consistent with the ISMF direction of the SI shell. The distribution of 
nearby interstellar Fe"*" with log A^(Fe"'")< 12.5 cm~^ is described equally well by 
a uniform distribution or an origin in spherical shell-like features. Higher column 
densities of Fe"*" (log A^(Fe"'")> 12.5 cm^^) tend to be better described by the 
pathlength of the sightline through the SI and S2 subshells. Column densities 
of the recombinant ion Ca+ are found to increase with the strength of the inter- 
stellar radiation field, rather than with star distance or total pathlength through 
the two magnetic subshells. The ion Ca^ can not be used to trace the distri- 
bution of local interstellar gas unless the spatial variations in the radiation field 
are included in the calculation of the ionization balance, in addition to possible 
abundance variations. The result is that a model of Loop I as composed of two 
spherical magnetic subshells remains a viable description of the distribution of 
nearby low density ISM, but is not yet proven. 

Subject headings: ISM: magnetic fields, bubbles, clouds, structure — cosmology: 
cosmic background radiation — Sun: heliosphere 
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Introduction 



The location of tlie Sun in the rim of the Loop I superbubble has been inferred from 
radio continuum data, kinematical data on the flow of local ISM away from the center of 
Loop I, data on gas-phase abundances in local ISM, and the coincidence of the velocity of 
ISM inside and outside of the heliosphere. Loop I is an evolved superbubble shell formed 
from stellar evolution in a subgroup of the Sco-Cen association, ~ 4 — 5 Myrs ago (e.g. 
de Geud Il992l : iFrischI Il995l . 1 19961 : iMaiz-ApellanizI 120011). Both t he original dimensions found 
for the Loop I bubble observed in 820 MHz (IBerkhuiisenl Il973l ) . and more recent studies of 
Heiled (jl998al Jbl. H98a,H98b) and IWoUeberu ( 120071 ). place the Sun in or adjacent to the rim 
of a magnetic superbubble shell for an assumed spherical geometry. The 1.4 GHz Wolleben 
study defines two magnetic subshells of Loop I, SI and S2, with magnetic pole directions 
differing by 90° ± 42°. Comparisons between the radio continuum filaments of Loop I and 
optical polarization data indicate that the radio filaments at distances of ~ 100 — 150 pc trace 
magnetic field lines, indicating that optical polarization is a suitable tracer of magnetic shells 
(H98a). Both the kinematics and abundance pattern of local interste llar material (LISM) 
suggest that the Loop I remnant has expanded to the solar location (IFrischI 1 198 11 1. LISM 
abundances of the refractory elements Mg, Fe, and Ca , show the charact eristic enhancement 
indicative of grain destruction in interstellar shocks ( iFrisch et al.lll999l ). Local interstellar 
gas, d < 50 pc, and dust flow away from the ce nter of Loop I at a best-fit velocity of ~ 18 
km s"^ in the local standard of rest (LSR, e.g. iFrisch et al.ll2009l ). The first spectrum of 
backscattered Lya emission from interstellar hydrogen inside of the heliosphere showed that 
the velocity of interste llar H° inside of the heliosphere is comparable to LISM velocities 
(lAdams fc FrischI 119771 ). Together these data suggest that the magnetic field and spatial 
configuration of the LISM can be us ed to test wheth er the Loop I magnetic superbubble has 
expanded to the solar location. The IWoUebenl (120071 ) model of the SI and S2 shells provides 
enough detail to make preliminary comparisons between LISM data and the properties of 
these shells. These comparisons provide interesting insights into the LISM properties, and 
support the possibility that local ISM within ~ 55 pc is dominated by the SI and S2 shells. 

Superbubble expansion into ambient ISM with equal magnetic and thermal pressures 
yields roughly spherical superbubbles during early expansions st ages when magnetic pres- 



sure is weak compa red to the ram pressure of the expanding gas (IMacLow fc McCraylll988 



Ferriere et al.lll99l[). and bubb les elongated along the ISMF during late stages of evolution 
(IHanayama &: Tomisakal l2006l ). The evolved shell is thicker near the ISMF equatorial re- 
gions, where field strengths are larger due to flux freezing, than the polar regions of the shell 
where thermal pressure provides the main support for the shell. In media where magnetic 
pressure is weak, e.g. the ratio of thermal to magnetic pressure (3 > 10, the evolved bubble 
is more symmetric. Supernovae in Sco-Cen Association subgroups have contributed to the 
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evolution of the Loop I superbubble during the past < 14 Myrs. The Loop I superbubble 
(and SI, S2) expanded in a medium with a density gradient, because the initial supernova 
occurred in the molecular regions of the parent Scorpius-Centaurus Association subgroups, 
while the subs equent bubble expansion occurred in the low density interior of the Local 
Bubble cavity (jFrischlll98ll . Il995l : iFuchs et al.ll2006l ). In this case the external plasma /3 may 
have varied irregularly across the expanding shell, so that the topology of the present day 
SI and S2 shells may deviate from axial symmetry as well as sphericity. 

The ISMF direction at the heliosphere provides the most direct measure of whether the 
Sun is embedded in the shell of the Loop I superbubble. Sev eral phenomena trace the field 



direction - the weak polarization of light from nearby stars (ITiiibergen 



hereafter F07), the flield direction in the SI subshell of Loop I ( IWoUebenl I2007D . the 3 kHz 



1982 



emissions from the outer heliosheath detected by the two Voyager satellites (iGurnett et al. 



Frisch 2007a 



20061 . F07), the observed angular offset between interstell a r H° and He° flow ing into the 



heliosphere ( iLallement et al.ll2005l : iPogorelov &: Zankll2006l : lOpher et al.l 120071 ). and the 10 
pc difference between t he distances of the solar wind termination shock detected by the two 
Voyager satellites (e.g. IStond 120081 ). The orientation of the plane midway between the hot 
and cold dipole moments of the cosmic microwave background is also within ~ 15° of the 
local ISMF direction (F07). Q 

This paper searches for evidence that the SI and S2 shells affect the distribution of 
nearby ISM within ~ 55 pc. The topology of the SI and S2 shells is discussed in ^ Section 
[3] shows that the direction of the ISMF at the Sun is consistent with the ISMF direction 
in the SI shell, similar to the location of the mid-plane between the cosmic microwave 
dipole moments, and consistent with the ISMF direction inferred from heliosphere models. 
The distribution of the ISM in the SI and S2 shells are compared to Fe"*" column densities 
towards nearby stars behind the shells (§!]). A similar comparison is made between the Ca"*" 
data and the SI and S2 shells, however Ca"*" column densities appear instead to trace the 
strength of the local far ultraviolet (UV) diffuse radiation field (§5]). An appendix outlines 
the ionization equilibrium of Ca"*". 



^More recently, IBEX has found that the ISMF interacting with the hehosphere forms a ribbon of energetic 
neutral atom emission as viewed from the Earth, and the ribbon traces regions where the ISMF direction 
within the outer heliosheath is perpendicular to the sightline (e.g. McComas et al. 20091 ). 
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Approximating the Three-Dimensional ISM Distribution in the SI and S2 

Shells 



WoUebenl (120071 ) has fit two separate spherical magnetic shells ('SI' and 'S2') to the low 
frequency (1.4 GHz and 23 GHz) polarized radio continuum, which must have a relatively 
local origin because of the dependence of Faraday rotation. The ISMF is assumed to be 
entrained in the expanding superbubble shell, with no deviation from spherical symmetry. 
The Sun is located in the rim of the SI shell, which is centered 78 ± 10 pc away at galactic 
coordinates i,b = 346° ± 5°, 3° ± 5°. The upwind direction of the fiow of local ISM past the 
Sun is within ~ 20° of of the SI shell center 3 The inner and outer radii of the SI shell are 
72±10 and 91±10 pc respectively. Wolleben described the SI magnetic field direction by two 
angles, the angle between the field direction and the NGP B^= —72° ±30°, and the rotation 
about the NGP Bg= 71° ± 30°. The S2 shell center is more distant (~ 95 pc) and centered 
at higher galactic latitudes {b ~ 37°) than the SI shell, with an ISMF direction near the 
north galactic pole. An alternate single-shell model for Loop I is also based on the H° shell 
and centers the feature at i,b= 320°, 5° (H98a). As a first approximation of shell structure, 
they are assumed spherically symmetric, although the ISMF is seen in filaments interacting 
with denser clouds for more distant regions of Loop I (H98a) A more detailed model of these 
evolved bubbles requires understanding the magnetic pressure. The parameters provided for 
the SI and S2 shells by Wolleben are detailed enough for comparison with observations of 
the LISM. 

A three-dimensional (3D) spherically symmetric model of the SI and S2 shells is created 
for comparison with the local interstellar magnetic field and distribution of interstellar Fe"*" 
and Ca"*". The 3D configuration is initially constructed in the frame of each shell such that 
the north pole is at the local zenith. The shell rims are filled with a uniform density of 
points. The shell model is then rotated to the galactic coordinate system by and B^, 
and translated to the shell center in the galactic coordinate system. The result is a model 
where a path through the shell measures the shell column density normalized to an arbitrary 
value, and the shell column density in any sightline varies according to the end point of the 
path. Fig. [H right, shows the ISM distribution in the SI and S2 shells for slices parallel to 
the galactic plane and for 10 pc-wide intervals of Z above and below the galactic plane. Fig. 
[H left, shows sections of the two shells at distances of ~ 30 pc using an aitoff projection. 
Negative Z- values are dominated by the SI shell (red), below Z~ 30 pc (latitudes b ~ 20°). 



^Comparisons between the LSR flow velocity and Loop I require using a somewhat uncertain velocity 
correction to obtain the LSR motion of the cloud. For the 'Standard' LSR, the local ISM flows at a velocity of 
-19.4 km s~^ from the direction of £,b=331°-5°, while an LSR correction based on Hipparcos star distances 
gives a bulk flow velocity of -17 km s^^, from ^,fe=2°,-5° |Frisch fc Slavinlbooel) . 
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At higher latitudes a given sighthne may sample either, or both, of the shells. The parameters 
of the 3D simulation place the Sun in the rim of the SI shell, but the uncertainties quoted 
by WoUeben also allow the Sun to be in the S2 shell. As additional data on nearby ISM 
become available, it should become possible to both constrain and test the S1/S2 models in 
more detail. 

The SI and S2 shells overlap in places. One example is the sightline towards the stars 
a Oph (HD 159561), located 14.3 pc aw ay at £,6= 35.9° , 22.6°, which has the strongest Ca+ 
line observed towards any nearby star (jCrawfordI l200ll ) . The SI and S2 shells coincide at 
a distance of 12 pc in this sightline, suggesting that the Ca"*" line towards a Oph samples 
a region where the S2 shell collided with the SI shell, possibly creating a shock so that 
recent grain destruction occ urred. Merging flows indu ce thermal instabilities that generate 
such filamentary structures (I Audit fc Hennebelldl2005l ). and the more distant H° gas in this 
sightline i s also filanientar y. A second example of possible interacting shells is the nearby Leo 
filament (ILauroeschl 120071 ) . The orange symbol in Fig. [1] (Z=25-35 pc) shows the location 
of a tiny cold (20 K) filamentar y (> 7° x 2° ) cloud in Leo, l ocated at £,h= 220°, 45° and 



at a distance of less than 42 pc (IMeyeii 120071 : ILauroeschl 120071 ) . If the cloud is at 40 pc, it 
is outside of both shells for the basic values for the shell distances and radii (i.e. without 
invoking any uncertainties). However if the cloud is nearby, or if the extreme values allowed 
by the uncertainties on the SI and S2 shells are invoked, this filament may form where the 
two shells collide. Merging flow s induce therma l instabilities that generate such filamentary 
structures as the Leo filament (ILauroeschl l2007l : I Audit &: Hennebellel l2005l ) , so the presence 
of this cold filament is consistent with the picture of the LISM as dominated by two shell 
features. 



3. The SI and S2 Shells and Local Interstellar Magnetic Field 



3.1. SI Shell and ISMF Direction 



Polarization by charged irregularly shaped interstellar grains yields polarization vectors 
that are parallel to t he ISMF, because the induced magnetic torques naturally align the grains 



(e.g. lLazarianll2000l ). There is a patch of dust towards the fourth galactic quadrant (£> 270°), 
mainly in the southern hemisphere and in the upwind direction of the interstellar gas flowing 
through the heliosphere, wh ere the ISMF wit hin 5-40 pc of the Sun has been traced by very 
weak optical polarizations ( iTinbergenl Il982l ). The Tinbergen data were acquired in the 
southern hemisphere during 1974, and northern hemisphere data during 1973 (J. Tinbergen, 
private communication). Tinbergen detected polarizations of > 0.017%, with la ~ 0.007%, 
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towards a few stars within 40 pc. Five of these stars § are close to the echptic plan e, and 



offset by up to A< 40° from the heliosphere nose towards positive ecliptic longitudes ( iFrisch 



20051 . l2007al ). The mean position angles for the three stars near the nose with P > 3cr, are 
PAg=33° ± 11° in galactic coordinates, and PAe=— 26° ± 11° in ecliptic coordinates (Table 
1). The position angle uncertainty is estimated by allowing Q and U to vary over ilex 
(Tabled]). The nearest star towards the nose is 36 Oph, at 6 pc, and it has a 2.5cr detection 
with P = 0.018%. The polarization position angle of 36 Oph (PAe=-19.9, PAg=39.5) does 
not differ significantly from the mean position angle of the three more distant stars with 3a 
polarizations in the nose regi on. The polar izations of the Tinbergen's sample, together with 



northern hemisphere data of iPiirolal ( 119771 ). are plotted in Fig. [2l Ar i alternative catalog of 



nearby star polarizations is the comprehensive catalog assembled by iLeroyi (119931 ): however 
it is mainly based on measurements with larger uncertainties, and is therefore less useful for 
identifying very low polarization levels. 

In order to test the relation between the ISMF in the SI shell and the Tinbergen 
polarization data, the ISMF direction for the SI shell was varied within the uncertainties on 
Bg and to find the direction that is the most consistent with the optical data. The best 
match to the optical polarization data was found for 71° and B^= —42°. The parts 
of the SI shell within 30 pc are compared to the optical polarization data in Fig. [2l for 
galactic (right) and ecliptic (left) coordinates, for stars within 50 pc; polarization vectors are 
plotted if P > 2.5cr. For the purpose of this figure, I use a center position for the SI shell 
of {i,b)= (351°, —2°) and radius 75 pc, so that the Sun is located 3 pc outside of the shell 
rim. The mean position angle of the SI magnetic field (for B^=-A2°) at the star locations 
is PAsi^E= —20 ± 5 in ecliptic coordinates, which is within the uncertainties of the stellar 
polarizations. A slightly larger value of exactly matches the mean position angles of the 
starlight polarizations, but violates the quoted uncertainties. Therefore, the ISMF directions 
from the Tinbergen data and SI shell configuration are consistent to within the uncertainties. 

The Tinbergen stars with the strongest polarizations are ~ 90° from the pole of the 
SI shell, which is consistent with the expectation of higher ISMF field strengths where the 
shell expansion is perpendicular to the ISMF direction. The best value for the SI magnetic 
field direction close to the Sun, derived from comparisons with these optical polarization 
data (§3]), is Bq= 71° and B^= —42°, corresponding to a local ISMF direction towards 
£,6= 71°, 48°. 



■^The five stars with strongest polarizations in the heliosphere nose region arc HD 161892, HD 177716, 
HD 181577, HD 155885, and HD 169916, with the first three stars showing polarization detections at the 3a 
level. 
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Heliospheric asymmetries are caused by interactions with the interstellar magnetic field. 
A widely used measure of the heliosphere distortion due to the ISMF, which is inclined 
by the angle a ~ 30° — 60° with respect t o the ISM flow vector, is the ob served offset 
between the inflowing He° and H° directions Jwittelbooi llallement et al.lboosi). The ISMF 



also s hifts the maximum Lya emission originating in the outer heliosheath (IBen-Jaffel et al. 



20001 ) ■ Correcting the He° and H° directions to a common observation epoch yields an 



offset angle separation of 4.9° ± 1° between the two directions. (The upwind direction 
of the He° flow in J2000 coordinates is A= 255.4° ± 0.5°, /3= 5.1° ± 0.2°, Witte, private 
communication). These directions deflne a position angle, which is PAe= —35° ±20° (PAg= 
26° ± 20°) in ecliptic (galactic) coordinates, respectively. Large uncertainties are quoted 
because the upwind direction of the H° flow through the heliosphere is not precisely deflned 
due to the ~ 50% of the interstellar H° lost to flltration in the hydrogen wall, the balance 
between radiation pressure and gravity affecting trajectories of H° atoms surviving to the 
helios phere interior, and the produc tion of secondary H° atoms inside of the heliosphere 

(e.g. 

£,6= 



Quemerais fc Izmodenovl |2002| ). For comparison, at the heliosphere nose location, 
3.5°, 15.2°, the SI shell with B(^=-42° gives a position angle PAe= — 15°. The position 
angle formed by the offset between He° and H° flowing through the heliosphere is marginally 
consistent with the SI shell direction at the heliosphere nose. The 10 AU difference in 
the termina tion shock distance found by Voyagers 1 and 2, in 2004 and 2007 respectively 
(IStonell2008l ). must be combined with the H°-He° offset to provide a more reliable constraint 
on models of the direction of the inters tellar magnetic fleld affecting the heliosphere (e.g. 
Pogorelov et"aDl2007l : lOpher et al.ll2007h . 



3.2. SI shell and the CMB Dipole Moment 



The great circle that is midway between the hot and cold poles of the cosmic microwave 
background (CMB) dipole pas ses within 5° of the interstellar He° upwind direction, and 
bifurcates the heliosphere nose (jFrischll2007al ). The point of closest approach to the nose is at 
£,6= 7.4°, 11.6°, where the position angle of the CMB dipole mid-plane in ecliptic coordinates 
is PAe= —11° ± 1° (uncertainties in the upwind He° direction, used to deflne the nose 



* Since this paper was originally submitted, a number of rec ent papers have appeare d that discuss the ISMF 
direction at the Sun, based on M HD heliosphere models ( e.g. iRatkiewicz et al. 2008: iPogorelov et al. 2008.') . 



IBEX data on the EN A Ribbon ( Schwadron et al. 2009: Funsten et al. 



2009f ). or both (jHeerikhuisen et al 



20101 ). There is an overlap between ISMF directions in these models and the uncertainties on Wolleben's 
ISMF direction for the SI shell. The overlap occurs for galactic longitudes in the range of 40° — 50° and 
galactic latitudes in the range of 23° — 42°. The models also predict an ISMF direction that is directed 
towards negative ecliptic latitudes. 
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position, are included in this uncertainty, Table 1). For comparison, the SI shell magnetic 
field direction at this location, for B^= —42°, corresponds to a position angle (ecliptic 
coordinates) of PAe= — 18°. This fact is mentioned here because the low-i multipole moments 



of the CMB show symmetries related to the ecliptic geometry (e.g. ICopi et al.ll2006[ ). so that 
the symmetry of the CMB dipole moment around the heliosphere nose should also be of 
interest. The upwind direction is 5° above the ecliptic plane, so any coincidence between the 
CMB multipole moments and the ecliptic geometry is tantamount to a coincidence with the 
heliosphere morphology and/or to the interstellar magnetic field that shapes the heliosphere. 
These ecliptic signatures on the CMB are not understood, but it is not unreasonable to 
postulate they arise fromj)rocesses related to the local interstellar magnetic field and its 
effect on the heliosphere. 



4. Comparisons between SI and S2 Shells and Distribution of Fe^ 



If the distribution of nearby ISM is determined by the SI and S2 magnetic superbubbles, 
then the SI and S2 shell morphologies should be imprinted on the strengths of interstellar 
absorption lines. Fig. [1] shows two views of the ISM associated with the SI and S2 shells. 
Interstellar absorption lines to wards stars within 55 p c show that most of the the LISM 
is warm, ~ 3,000 — 12,000 K (IRedfield fc LinskyI 120041 ) with low average spatial densities, 
< 0.1 cm~^. Exc eptions are the tiny dense clouds occasionally seen in Na° and H° absorp- 



tion ( iMeyerl 120071 ). Models of the radiative transfer properties of the circumheliospheric ISM 

0.07 cm~^, and with 



,-3 



show a partially ionized, low density cloud, n(H°)~ 0.20 cm 
temperature ~ 6, 300 K de termined from interstellar He° inside of the heliosphere ( IWitte 



2004J : ISlavin fc Frischl l2008l SF08). Therefore, a suitable tracer of the SI and S2 shell mor- 
phologies should be abundant, insensitive to cloud ionization, and undepleted. There are 
no available data sets that meet all three requirements, therefore the criteria that the el- 
ement be undepleted is dropped. The best ele ment for this study is then Fe"*", which has 
been measured towards ~ 27 stars within 56 pc (ILehner et al.ll2003t IRedfield fc Linskyll2002 



Kruk et al.ll2002l ). Iron is predominantly singly ionized in the cloud around the heliosphere, 
with neutral and Fe"*"^ together containing less than 3% of the Fe atoms (SF08). A more 
difficult aspect of using Fe"*" to trace absolute ISM densities is the factor of ~ 6 — 40 difference 
in the gas-phase Fe+ abundances between dense cold and warm tenuous clouds, due to dus t 
grain destruction by interstellar shocks including in local regions (e.g. ISlavin et al.l 120041 ). 



^In results published after the submission of this paper, it has been shown that the ISMF interacting 
with the hel iosphere controls the flow of nanometer-sized interstellar dust grains around and through the 
heliosphere ( Slavin et al. 20091 ). 



9 



The alternative common element that traces both neutral and ionized gas is Mg"*", however 
it has similar abundance variations as Fe+ and the Mg"*" h and k lines may be more saturated 
than the Fe+ lines. Therefore, I use Fe"*" column densities to trace the distribution of local 
ISM. Three nearby stars are omitted from this discussion because they have known debris 



disks (HD 215789, HD 209952, HD 216956. ISu et al.ll2006[ l. 



Comparisons between log A^(Fe^) and the column density of shell gas (SI and S2) in 
front of each star, versus a comparison between iV(Fe"'") and the star distance, provide useful 
insights (Fig. [3]). The shell column density towards each star represents the sum of the 
densities through the parts of the SI and S2 shells foreground to the star, normalized to 
an arbitrary value. This pathlength was constructed by assuming a column width of ±5 pc 



in or der to smooth out uncertainties in the intrinsic shell parameters (Table 1 in IWoUeben 



20071 ). Because of this smoothing, stars within 6 pc are omitted from Fig. El The bar 
at the bottom of the figures shows the column density range for stars within 6 pc. For 
low column densities, log A^(Fe+)< 12.5 cm~^, Fe"*" column densities tend to increase with 
star distance. For higher column density sightlines, log A^(Fe^)> 12.5 cm~^, A^(Fe"'") clusters 
more tightly around the shell pathlength (e.g. column density) than around the star distance. 
The exceptions to both comparisons are the stars HD 120315 (a low column density, high- 
latitude star that should sample a long pathlength through the shells) and HD 80007 (a high 
column density low-latitude star, with a path that is tangential to the SI shell). The Fe"*" 
components in a sightline are summed together for this comparison. The dashed lines in 
Fig. Oshow a linear fit between the column densities and the ordinate. The individual stars 
are listed in Fig. HI where the relative column densities of the SI and S2 shells towards each 
star are shown. 

A second property of the distribution of Fe"^ in these figures is that stars with galactic 
longitudes > 180° tend to have larger Fe"*" column densities, at a given distance, than stars 
in the opposite hemisphere. This effect is not seen in D° or H°. One possible reason is that 
the Fe abundances differ between the two hemispheres. The problem with this explanation 
is that it requires abundance variations over spatial scales of several parsecs, in relatively low 
velocity ISM (< 20 km s~^ LSR), with the variations ordered by the arbitrary coordinate of 
galactic longitude. An alternative explanation is that the Fe"*" lines for 180° include cool 
unresolved clumps of ISM. The ISRF towards l> 180° is significantly larger than towards 
180° (^|5]), so that ISM for £< 180° experiences reduced heating because of the absorption 
of H-ionizing photons, which may allow unresolved ISM clumps to coexist with warmer gas 
at the same velocity. Once the ionization of local ISM is better understood, the morphology 
of the SI and S2 shells can be adjusted to better represent the actual distribution of the 
local ISM. 
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The stars within 6 pc show total column densities log A^(Fe"'")= 12.14 — 12.88 cm~^. Of 
these stars, the strongest lines are towards HD 155885, HD 165341, and HD 187642 in the 
upwind direction. 

The distributions of SI and S2 shell material towards each star are shown in Fig. |H 
Stars with high (low) column densities are plotted as large (small) symbols for each direction. 
The most evident property is that the SI shell dominates the ISM towards the southerly 
stars. These can be used to predict whether a star (or exoplanet system for example) 
will be embedded in a cloud-like feature or in the Local Bubble plasma. The assumed 
spherical morphology for SI and S2 leads to the predictions that the white dwarf HD149499B 
(WD1634-573) is embedded in in the relatively denser gas of the shell, while the white dwarf 
WD1615-154 will be embedded in the low density Local Bubble plasma. 



5. Ca+ and Electron Densities in the SI, S2 Shells 



Interstellar Ca^ is a recombinant species that traces the electron density as well as 
abundance variations. Because Ca+ is formed through recombination, it is a proxy for 
the electron density in nearby low density gas providing that abundance variations and the 
radiation field are understood. If the ISRF and electron densities are uniform throughout the 
SI and S2 shells, then Ca"^ column densities would show a similar dependence on pathlength 
through the shells as seen for Fe"*". Interstellar Ca"*" column densities are plotted against the 
pathlength through the SI and S2 shells (Fig. [5l left) a nd star distance (Fig. [5], right), using 
data from (IFrisch et al.ll2008l . l2002t IWelty et al.l 119961 ). The Ca"*" column densities do not 
correlate with either the star distance or pathlength through the SI and S2 shells. Sightlines 
where the Ca"*" column density is an upper limit are not included in this comparison. 

Higher Ca+ column densities are found for stars with l> 180°, as was seen for Fe"*". 
The mean Ca^ column density is ~ 40% higher for stars in Quadrants III and IV, £> 180°, 
compared to stars with £< 180°. The difference becomes a factor of two if the anomalously 
strong Ca+ line towards a Oph is ignored. For 180°, <A^(Ca+)>= 3.7 x 10^° cm~^ (23 



stars). Omitting a Oph, which has the strongest known Ca+ for nearby stars (e.g. I Crawford 
200lh . gives a mean for the ^< 180° sample of <A^(Ca+)>= 2.3 x 10^°. For ^> 180°, 



<A^(Ca^)>= 5.1 X 10^° cm ^ (23 stars). No s imilar effect is see n in H° (or D°) column 



densities towards nearby stars (based on data in lWood et al.ll2005l ). For stars within 50 pc, 
the mean H° column density does not vary between the £< 180° hemisphere and the £> 180° 
hemisphere, and both show a mean value of A^(II°)~' ^ ^ ^ ^"^^^ 
sizes of 25 and 27 stars. 



1.5 X 10 cm for respective sample 



11 



A different picture emerges when Ca+ (and Fe+) column densities are compared to 
the far UV radiation flux at the star. The highest diffuse far UV fluxes are seen towards 
stars in the third and fourth galactic quadrants, £> 180°, because of the low ISM opacity 
i n the Local Bubble interior, hot stars in the Scorpius-Centaurus Association, and a Vir 



(IGondhalekar et al.lll980l : lOpal fc Welledll984 Go80,OW84). Fluxes at 975 A are given for 



the 25 brightest stars, based on a survey by the STP 72 — 1 satellite in the 910-1050 A band 
and flux models (OW84). At least 97% of the local flux at 975 A is provided by the stars 
with i> 180°. Using the 975 A data for the 25 brightest stars (OW84), the 975 A flux was 
calculated at each star and is plotted against Ca+ column densities (Fig. [6l left), and Fe"*" 
column densities (Fig. |6l right). The Fe^ and Ca+ samples are different, with more than 
half of the Fe+ stars within 20 pc, while most of the Ca"*" stars are beyond 20 pc. All of 
the stars in the Fe"*" data set with 975 A flux levels larger than 6.8 x 10^ photons cm~^ s~^ 
A~^ are in galactic Quadrants III and IV, and these stars tend to have larger Fe"*" column 
densities. The ISRF gradient would affect iV(Fe^) only through larger column densities of 
HII, since Fe^ dominates in both ionized and neutral diffuse ISM. Nearby H II gas in regions 
with t> 180° would explain Fe"*" variations, without violating the H° constraints (which are 
determined from a similar star sample). 

The dependence of Ca+ column densities on the 975 flux is less simple because Ca+ is a 
trace species formed by recombination, with an ionization potential of 11.87 eV versus 13.60 
eV for H°. In the cloud around the Sun, n ~ 0.2 cm'^ and Ca++/Ca+=63 (SF08). The ISRF 
gradient near the Sun (Fig. E]) affects A^(Ca+) several ways. Higher radiation fluxes lead to 
higher electron densities, increasing Ca++^Ca+ recombination, and the overall H"*" fraction 
would increase. Higher radiation fluxes also increase the Ca"*" photoionization rate, but this 
effect does not appear to be dominant. In the local ISM, photoionization appears to dominate 
over collisional ionization. Radiative transfer models of the ISM surrounding the Sun match 
available ionization data such as the Mg+/Mg° ratio (SF08). Low ob served Ar° abundance s 



towards nearby stars also indicate the dominance of photoionization ( ISofia fc Jenkinslll998l ). 
The radiation flux at 1044 A capable of ionizing Ca+ is traced by the 975 A radiation field, to 
within ~ 10% (Go80,OW84), so that Ca+ ionization rates should increase with the 975 flux. 
The increase of A^(Ca+) with radiation flux is predicted by the photoionization equilibrium 
of Ca+ (see the appendix). Predicted A^(Ca~^) values are plotted against the ISRF for three 
different total H column densities, A^(H)=A^(H°)+A^(H"'"), in Fig. |6l The observed increase 
of A^(Ca"'") with higher fluxes is consistent with Ca^ photoionization, and indicates that 
higher electron densities and H II column densities are both significant factors in the Ca"*" 
line strengths. 
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Discussion and Conclusions 



The discussions in ttiis paper are based on a searcti for evidence of tlie SI and S2 stiells in 
local ISM data. The SI and S2 shells are assumed to be spherical and complete. Such simple 
assumptions are justified only in the initial stage of probing the ISM distribution associated 
with a superb ubble shell that has column densities too low for H° 21-cm measurements, 
and that has evolved into a very low density region of space. Several studies model the 
formation of the Local Bubble in terms of the energy iniected into the ISM by su pernovae 
in the Sco-Centaurus Association (jde Geuslll992l : lFrischlll998l : lMaiz-Apellanizll200ll ). but the 
connection between the Loop I radio emission and very local ISM has never been established. 
The SI and S2 shell models provide a basis for testing this connection. 



Early optical polarization data ( lTinbergerull982l ) indicate that the ISMF direction close 
to the Sun agrees with the SI shell ISMF direction once the uncertainties in WoUeben 's 
angle (2007) are included. In principle Fe"*" can be used to trace the ISM distribution, since 
is arises in both neutral and ionized gas. Detailed comparisons between Fe"*" line strengths 
towards nearby stars and the projected pathlength through the SI and S2 shells towards that 
star support, but do not prove, that the lines arise in shells. Both the Fe"*" and Ca"*" data 
indicate that the portions of the SI and S2 shells with £> 180° will be more highly ionized 
than in the opposite hemisphere. The data are not sufficient to distinguish an ionization 
gradient from an abundance gradient. For this reason, the shells are better traced using 
ions with first ionization potentials less than 13.7 eV. Heating by Lya radiation accounts for 
~ 66% of the heating of the circumheliospheric ISM, so shell regions exposed to the highest 
radiation fiux should also be warmer (an effect not explicitely included in the Ca equilibrium 
discussion in the appendix). 

These models of the SI and S2 shells assume spherically symmetric forms, which may be 
viable only for low density sightlines where magnetic and thermal pressures are comparable. 
Interstellar data are compared mainly to t he SI she ll , whi ch has the most favorable geom- 
etry for surrounding the Sun according to IWoUebenl (120071 ). The Tinbergen data suggest a 
slow increase in polarizations with distance, and the position angles t owards the no se are 
consistent with the optical polarization of more distant stars in Loop I (jFrischll2007al ). 



The most distant shell regions, in the galactic center hemisphere, have expanded into 
the high-extinction g as beyond ^ 100 pc that is ass ociated with t he Sco-Cen Association 
(see e.g. Figs. 1,3 in I Chen et al.l ( 119981 ) or Fig. 2 in iFrischI ( l2007bl )) and show pronounced 
magnetic filaments (H98a) rather than a spherical shell geometry. Heiles points out that 
the synchrotron-emission ridges of Loop I follow the distortion of the nearby global ISMF, 
as traced by polarization data, and that Loop I is not a shell for the high density regions. 
The global ISMF within several hundered parsecs is directed towards £~ 80° — 88°. The 
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expansion of the nearest portions of the SI shell in a uniform field would yield an ISMF 
close to the Sun directed upwards with respect to the galactic plane. 



Hanayama &: Tomisakal (120061 ) model the properties of a magnetic superbubble ~ 3 Myrs 



old and for the strong ISMF case, where magnetic pressure dominates thermal pressure by 
a factor of ~ 2. The superbubble cavity is elongated in the direction parallel to the ISMF, 
where the shell is thinner. The shell is thicker and more extended in the radial direction, 
~ 90° from the magnetic pole. The region of strongest polarization for the Tinbergen sample 
is ~ 90° from the magnetic pole. Non-uniform shell expansion, or a 'wrinkled' shell, could 
explain the anomalous sightlines towards HD 120315 and HP 80007 (se e Fig. |3] compared 



to Fig. H]), while the tiny dense nearby Leo cloud ( iMeyer et al.l l2006l . Lauroesch private 



communication), and the exceptionally strong Ca"*" line towards Rasalhague (HD 159561, a 
Oph) may indicate a region of merging shells. 

The column densities for Fe"*" and Ca"*" are generally weaker for sightlines with i< 180°, 
and stronger for stars with £> 180°. This effect may be either from the distribution of 
ionized gas, or abundance variations for Fe and Ca. The effect is seen over small spatial 
scales of ±10 pc. If the variation is due to abundance differences, then the ISM close to the 
Sun would have two different histories, although the flow velocities are similar. 

The conclusions of this comparison between the SI and S2 shells with LISM markers 
can be briefly summarized: 



The IWollebenI (120071 ) description of the Loop I polarized radio continuum in terms 
of two shells, SI and S2, is viable and has sufficient detail to be tested against ob- 
servational data. For example, when the the SI and S2 shell paramete rs are varied 
within the allow ed uncertainty range, the nearby cold gas filaments in Leo (jMeyerll2007l : 
LauroeschI 120071 ) are seen to be produced where the two shells merge or collide (f|2]). 



The SI shell magnetic field direction of IWoUebenl (|2C)07), with Bn',=- 42°, matchs the 
ISMF direction derived from older polarization data ( lTinbergerull982l ) of nearby stars 
near the ecliptic plane and heliosphere nose, but offset by up to A ~ +40° from the 
heliosphere nose. The ISMF direction implied by the SI shell and polarization position 
angles together is directed towards £,6= 71°, 48°. 

For low column densities, log A^(Fe~^)< 12.5 cm^^, the strength of the A^(Fe"'') is 
better described by the star distance. For higher column densities, log A^(Fe"^)> 12.5 
cm~^, the strength of the A^(Fe"'") is better described by the pathlength of the sightline 
through the SI and S2 shells. This result is based on a limited number of stars (< 25) 
and requires confirmation using a larger data set (§!]). 
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• The illumination of the SI shell by the strong diffuse far ultraviolet interstellar radiation 
field in Quadrants III and IV, l> 180°, explains the higher column densities observed 
for Fe+ and Ca"*" in these galactic quadrants (§11 An appendix evaluates the 
ionization equilibrium of Ca"*" features spaced around the shell, and shows that the 
local radiation field strength regulates the Ca^ absorption line strengths. 

• The ISMF direction at the heliosphere nose is within 10° of the angle of the great circle 
that is midway between the hot and cold hemispheres of the CMB dipole moment, and 
that also bifurcates the heliosphere nose. 

• The S1/S2 shell model can be used to predict whether a star, or exoplanet system for 
example, is embedded in a cloud or in the Local Bubble plasma (S]). The reverse is also 
true, that measurements of astrosphere properties will help constrain the distribution 
of ISM associated with the shells. 

• This scenario describing the influence of the magnetic superb ubble SI and S2 shells 
on the local ISM, and as the origin of the interstellar magnetic field at the Sun, is 
consistent with available data, but does not yet prove the S1/S2 model. Two kinds 
of data are required to substantiate this picture: (1) Additional UV observations of 
tracers of both neutral and ionized interstellar gas, e.g. Fe+, Mg"*", and Mg° features. 
(2) Measurements of nearby weak interstellar polarizations at 0.01% levels or better. 



The author would like to think NASA for research funding, in the form of grants NAG5- 
13107 and NNG05GD36G to the University of Chicago. 



A. Ca+ Equilibrium 

The Ca"^ column density can be determined from the assumption of photoionization 
equilibrium between Ca"*" and Ca^"*", for the radiation flux level F\ at each star. The highest 
fluxes of diffuse far UV radiation are seen towards the third and fourth galactic quadrants, 
^> 180°. Some self- shielding of the ISM in the two shells may occur, but the very low opacity 
of the Local Bubble interior suggests that the radial dependence of the ISRF from the 
25 brightest far UV stars is more germane for understanding local ISM ionization, and the 
recombinant species Ca+ that tracks the ionization. For the te mperature rang e considered 



here, 3,000 — 15,000 K, collisional ionization is insignificant (jPottaschl Il972[ ). The Ca"*" 



equilibrium depends on the Ca+ photoionization rate r23, the recombination rate from Ca"^ 
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to Ca"^ 0^32, the electron density Ue, the Ca abundance Aca, and the total hydrogen density 
(H° + H+): 



Ca^ = 032 ne Ca^ 



(Al) 



or 



123 Ca+ = «32 He (^Ca (H° + H+) - Ca^ 



(A2) 



We define the ratio of the ionization and recombination rates at each star as 



*^'Mod26 * 



F2& 

, 2^-0.8 



T- 

_\J_-j0.8 
^26 



0.8 



(A3) 



for radiation flux F\, electron temperature Tj, 032 ~ j^-o s ( jShuU fc van Steenberg]|l985l ). 
and after parameterizing $ in terms of the properties of the circumheliospheric ISM (CHISM) 
at the solar location based on Model 26 in SF08. For Model 26 in SF08, ne=0.0654 cm~^ (and 
n(H+)=0.0554 cm'^), Ca++/Ca+=63.5, Tchism = T26=6320 K, giving $Mod26 = 4.15 using 
eq. IA2[ The ionization edge of Ca"*' is at 1044 A. For the ratio Fi/F2Q, F[ is approximated 
by the total flux at each star i from the combined distance-corrected radiation flelds of 
the 25 brightest stars at 975 A (base d on fluxes as measured at the Sun, with no opacity 
corrections, fromlOoal fc Wellerlll984j) . The normalization factor F2Q is the 1044 A flux from 

80, 000 photons cm"^ s"^ 



Gondhalekar et al. 



(119801 ) as used for Model 26 in SF08, or F26 
A~^. With this scaling, 6 = 4.74 x 10~^ in eq. I A3 1 Calcium abund ances appear to vary by 
a factor of ~ 40 between cold and warm clouds ( IWelty et al.l 119991 ). I use the typical warm 
diffuse cloud calcium abundance of Aca. = 2.2 x 10~^ calcium atoms per hydrogen atom. 

iV(Ca"^) then becomes: 



iV(Ca^ 



Ue A 



Ca 



tot) 



b Fi T"-« + n, 



(A4) 



The observed relation between Ca"*" column densities and the 975 A flux seen in Fig. E] 
is compared to predicted values of A^(Ca"'') determined from eq. IA41 for cloud temperatures 
in the range of 2,500 and 15,000 K and electron densities in the range of ne=0.01-0.15 cm~^ 



and positiye detections of Ca "*" for stars within 55 pc (using Ca"*" data from iFrisch et al. 



2002. 


2008; 


Weltv et al. 


1996) 



^The Ca"*" column densities represent the sum of all components towards each star. The stars are HD 
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temperature and electron density, Tj in eq. A3, and for that I use the somewhat arbitrary 
relation T = 200 n{e)~^'^^ in order to account for the increased cooli ng resulting from 
the c ollisional excitation of C"*" and 0° fine-structure levels by electrons (lYork &: Kinahan 



19791 ). This assumption is required to fully specify the ionization equilibrium. The lines 



labeled 20.00, 19.5, and 19.0 in Fig. [6] show the predicted Ca"*" column densities for these 
assumptions and log A^(H)= log A^(H°+H+) = 20.00, 19.5, and 19.0 cm'^. If the gas is 
clumpy, or conditions differ substantially from the CHISM gas, these estimates will break 
down. In the absence of a full 3D model of opacity over several hundred parsescs, more 
detailed comparisons between A^(Ca+) and the ISRF require additional data on the cloud 
temperature, H density, or ionization. However this comparison illustrates that reasonable 
assumptions for the parameters required to calculate the ionization equilibrium of Ca"*" yield 
reasonable predictions for the sensitivity of A^(Ca"^) to the far UV radiation field. 
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Table 1. Position Angles 



Item 


Galactic 


Ecliptic 




Coords. 


Coords. 




PAG(deg) 


PAE(deg) 


Position angles for 3 upwind stars with P > 3a^^^ 


33 ±11 


-26 ± 11 


Mean SI shell i?-field towards 3 upwind stars*^^-* 


39 ±6 


-20 ±5 


H°- He° offset at heliosphere nose {i,b= 3.5°, 15.2°) 


26 ±20 


-35 ± 20 


SI shell 5— field orientation at heliosphere nose 


46 


-15 


Direction of CMB dipole midplane at (7.4°, 11.6°; 


50 ±1 


-11 ± 1 


SI shell 5-field orientation at (£,6)~ (7.4°, 11. 6°)^^) 


43 


-18 



(^)These th ree stars are HP 1 61892, HD 177716, HD 181577, and the polarization 
data are from iTinbergenI (1l982l ). The nominal uncertainties on the position angles 
are obtained by letting Q and U vary over la measurment uncertainties of ±0.007%. 



{s)Xhe parameters for the SI shell are given in IWoUebenl ( 20071 ) 



C^^The H° inflow directionjs giyen in iLallement et all f l2005f ). The He° B1950 
inflow direction is given in lWittd (120041 ) . and must be corrected to J2000 coordinates 

(^^This angle is calculated for Be= 71° and -42°, an d is within the ±42° 



uncertainties on the SI shell ISMF direction ( IWoUebenl 120071 . after combining the 
±30° uncertainties on each angle in quadrature.) 

(£^)'Xhis location is the nearest point of the CMB dipole mid-plane to the helio- 
sphere nose direction, as defined by the interstellar He° flow through the heliosphere. 
Uncertainties in the He° fiow direction that defines the heliosphere nose are also in- 
cluded in thes e comp arison uncertainties. The CMB dipole directions are given in 
Bennett et aP (|l996h . 
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Fig. 1. — Left: The parts of the SI and S2 shells about 30 pc from the Sun. Right: 
Distribution of points in the SI (red) and S2 (blue) shells for different distance intervals 
Z (pc) above and below the galactic plane. The galactic center is at y=0 and x=infinity, 
beyond the plot right. The direction of galactic rotation is at x=0, y=infinity, beyond 
the plot bottom. The individual points trace the range of ISM distributions in the shells 
calculated with the assumption that each shell is spherically symmetric, and that the ISM is 
uniformly distributed. The dotted gray lines show the distance corresponding to 55 pc from 
the Sun for the midpoint of the Z-interval. The SI shell distribution is for a shell centered 
at i,b= 346°, 3° and 78 pc, with a radius of 82 pc, and with a total rim thickness of 19 pc. 
The S2 shell distribution is for a shell centered at i,b= 347°, 37° and a distance of 95 pc, 
with a radius of 75 pc and total rim thic kness of 24 pc. The orange st ar shows the location 
of the nearby cold filament towards Leo (IMeyerl 120071 : iLauroeschl |2007| ) . 
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Fig. 2. — The magnetic field associated with the parts of the SI shell within 30 pc is 
plotted in ecliptic coordinates (left) and galactic coordinates (right) for an aitoff projection. 
The ecliptic plot is centered at A = 2 70°, while the ga lactic plot is centered at £=0°. The 
parameters for the SI shell given in IWollebenI (120071 ) have been varied within the range 
of allowed uncertainties to yield the best match to the iTinbergeru ( 119821 ) polarization data 
towards stars in the heliosphere nose region. The dark and light blue dots show the inflow 
directions of interstellar H° and He° into the heliosphere. The SI subshell parameters used 
in the above figures correspond to a shell center at (^,6) = (351°, —2°) and 78 pc, shell radius 
of 75 pc, and magnetic field angles B0=71°, B^=-42°. The dots show stars within 50 pc 
with polarization data, and the red bars sho w polarization vectors for stars where P^> 2.5a 
(ITinbergeru Il982l : iPiirolal 119771 : iFrischI l2007al ) . The longitude increases towards the right for 
both figures. 
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Fig. 3. — Column densities of Fe"*" are plotted against the total predicted pathlength through 
the SI and S2 shells foreground to the stars (or equivalently column densities, in arbitrary 
units) through the SI and S2 shells (left, see ^ and against the star distance (right). 
The symbols indicate stars with stars 6-56 pc away and with galactic longitude £> 180° 
(open squares) and i< 180° (filled dots). The horizontal bar at the bottom shows the 
column density range for the stars within 6 pc that are not plotted. The vertical line 
separates stars with log A^(Fe"'')< 12.5 cm~^ versus > 12.5 cm~^ where column densities 
tend to increase with distance or pathlength through the shells, respectively. The stars are 
listed in Fig. |H The dashed line shows a first order fit between the total column densities 
and pathlength (left), and total column densities and star distance (right). The individual 
velocity components for each star are added to give the total iV(Fe"'") plotted on the abscissa. 
The highest Fe"*" column density in this sample is towards /3 Car (HD 80007), 34 pc from the 
Sun in the direction i,b= 286°, —14°. According to Fig. 1, right, and Fig. 4, left, this star 
samples a tangential pathlength through the SI shell. 
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Fig. 4. — The distributions of ISM in the SI shell (red) and S2 shell (blue) are shown 
towards each star in Fig. |31 based on a spherical morphology for each shell (see text). The 
distributions are plotted for arbitrary density units. The sum of the SI and S2 densities is 
plotted as the thick gray dotted line. The turquoise symbol shows the distance of the star 
(labeled by name and coordinates), and small (large) symbol sizes show log A^(Fe"'")< 12.5 
cm~^ (> 12.5 cm~^) towards the star. 
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Fig. 5. — Left: Column densities of Ca+ are plotted against the total predicted pathlength 
(in arbitrary units) through the SI and S2 shells sampled by each star. The stars are 5-55 pc 
away and are coded by whether the galactic longitude is £> 180° (open squares) or £< 180° 
(filled dots). The dotted line shows a first order fit between the total column densities and 
pathlength (left), and total column densities and star distance (right). The first two galactic 
quadrants, i= 0° — 180° show generally smaller Ca"*" column densities than do stars with 
i= 180° — 360°. The individual velocity components for each star are added to give the 
total A^(Ca+) plotted on the absicissa. Right: The same set of stars are plotted against the 
distance of the star from the Sun. 
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Fig. 6. — The interstellar radiation flux at 975 A at the location of each star is compared 
to Ca ^ column densities (left) and Fe"*" column densities ( right). The Ca" ^ data are drawn 
from JFrisch et al.lbonsi 120021 ). and the ISRF data is from lOpal fc Wellej jl984| ). Symbols 
indicate whether the galactic longitude of the star is less than 180° (filled), or greater than 
180° (open). The dotted curves show the predicted Ca"^ column densities for the labeled 
total hydrogen column densities (H°+H"'") for the assumption of photoionization equilibrium, 
and assumptions for the Ca abundance (2.2 x 10~^ per hydrogen atom) and temperature 
sensitivity to Ue (see appendix). The two stars with log A^(Ca"'")~ 10.8 cm~^ and low flux 
levels of ~ 4.8 x 10^ photons"^ cm~^ s^^ A^^ are between 50 and 56 pc away. 



